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Tangential and radial angles of muons in EAS, a useful concept in investigation of the muon production height, 
can be used also for the investigation of the muon momenta. A parameter £, being a combination of tangential 
and radial angles, is introduced and its possible applications in investigation of muons in showers are presented. 


1. Introduction 

EAS experiments, combining field detector ar¬ 
rays - being capable in measuring the shower di¬ 
rection - with muon tracking detectors, can be 
used for determination of the mean muon produc¬ 
tion height (MPH), a parameter sensitive to the 
primary mass. Muon tracks in air shower are, in 
general, not coplanar with the shower axis, which 
have led to the concept of tangential and radial 
angles Q]. However, as it will be shown in the 
following, these two independent observables may 
provide a possibility to investigate the muon mo¬ 
mentum components in showers. A good exam¬ 
ple of the EAS experiment, where the ideas pre¬ 
sented in this paper can be used is KASCADE [|J 
with its recently introduced large Muon Tracking 
Detector ||. We are using data simulated with 
CORSIKA |Q ver.5.644 and QGSJet hadronic in¬ 
teraction model for all discussions and examples 
given below. 

2. Tangential and radial angles 

As it is shown in Fig. [j], shower direction and 
the location of the muon detector define two per¬ 
pendicular planes: tangential and radial ones. 

Tangential angle r and radial angle p are angles 
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Figure 1. Definitions of tangential and radial an¬ 
gles. 


between the shower direction and the orthogonal 
projections of the muon track onto the tangential 
and radial planes, respectively. 

The measured value of r reflects the amount of 
muon scattering in the atmosphere and any rel¬ 
evant absorber/detector material, together with 
a possible displacement of the muon production 
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place from the shower axis. Its distribution (Fig. 
|^) is symmetrical around zero and gets narrower 
with the increase of the muon momentum, what 
one could expect from the quantity related to the 
multiple scattering. 



Figure 2. Dependence of the tangential angle on 
muon momentum. 
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Figure 3. Radial angle distribution. 


A typical distribution of the radial angle is 
shown in Fig. [| Its value depends on the muon 
parent’s (meson’s) transverse momentum and the 
distance of the muon to the shower axis at the ob¬ 
servation level. In real experimental conditions 
muons are sampled within certain range of dis¬ 
tances from the core. The shaded distribution is 
shown for a typical distance 20 m - 180 m from 
the core. The large p values are removed. 

3. Parameter £ 

As it is seen from Fig. |l], using simple 
trigonometry and substituting tangent for an an¬ 


gle, what is correct within 5% error up to 0.4rad, 
one can write the following expressions for r and 

P- 

T 9* £2 = £* x sin IPOB 

PD p|| 

P 2* £4 = a x cos LPOB 

r PD p|| 

Let us define the new parameter £ as follows: 

c = Vt 2 + p 2 = §■ ; 

which is valid for r < 0.4 rad, and p < 0.4 rad. 

So, measuring in an experiment angles of a 
shower and of muons and calculating p and r, 
by means of £ one gets the possibility to investi¬ 
gate the muon momentum space. The precision 
of C depends on the precision of angular measure¬ 
ments in the air shower, and, in particular, on the 
precision of muon tracking. It is known, that the 
tracking accuracy improves when the high energy 
muons only are considered. Therefore, one would 
like to have as large content of high energy muons 
in the analyzed data sample as possible. 
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Figure 4. Muon momentum distribution in show¬ 
ers initiated by proton primaries with the energy 
Uo=10 15 eV. The shaded area shows the effect of 
the experimental threshold lGeV/c and the lim¬ 
ited distance to the shower core: 20 - 180 nr. 


A typical muon momentum distribution in 
showers initiated by proton primaries with 
£b=10 15 eV is shown in Fig. [|. In all sim¬ 
ulated showers (solid line) « 30% are muons 
with momentum larger than 10 GeV/c. In real 
experimental conditions, where a threshold, eg. 
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1 GeV/c, is present and muons are registered 
in a restricted range of distances to the shower 
core, one gets the shaded distribution in Fig. |[ 
In this sample, with significantly reduced statis¬ 
tics, there is already « 45 % of muons with p^ 
> 10 GeV/c. However, one would like to increase 
high energy muon content in the data sample even 
more, and for this purpose the £ parameter can 
be used. 

Indeed, as seen from Fig. ||, £ shows strong de¬ 
pendence on the muon momentum. While at low 
momenta the full spread of £ values is observed, 
its distribution becomes more confined to the low 
values with the increase of the muon momentum. 



<£> 



Figure 5. Dependence of £ on the muon momen¬ 
tum. In the lower panel a mean value of £ is 
plotted. 


The mean value of £ decreases with momentum 
exponentially, what is shown in the lower panel of 
Fig. ||. So, it is possible to use in the analysis the 
cut on the £ value and, in this way, to increase the 
high energy muon content in the analyzed sample. 

In Fig. the ratio of muons with momentum 
larger than 10 GeV/c to all muons in a sample 


(above the threshold of 1 GeV/c) as a function 
of the applied cut on £ value is plotted for two 
primary proton energies. Only muons in the dis¬ 
tance between 20 and 180 m to the shower core 
were considered. It is seen that choosing for the 
analysis muons e.g. with £ < 0.02 results in sig¬ 
nificant increase of the high energy muon content 
in the data sample. 
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Figure 6. Relative content of muons with mo¬ 
mentum larger than 10 GeV/c in the whole data 
sample (muons above 1 GeV/c) as a function of 
applied £ cut. 


Another application of the £ parameter is 
rather straightforward. Its value is nothing else 
but a tangent of the muon angle versus the shower 
axis. So, knowing the distance to the shower 
core a simple triangulation allows to reproduce 
the muon production height. It is, off course, 
assumed that muons stem from the shower axis, 
what is generally not true. But, it turns out, that 
such approximation works well. 

In Fig. [i] the true muon production height dis¬ 
tributions (zq) are compared with the calculated 
ones (MPH), using the just mentionned proce¬ 
dure, for two primary species, H and Fe with en¬ 
ergy £’o=10 15 eV. In these calculations a cut on 
£ < 0.02 was applied. One sees a good agreement 
of both curves. The mean production height val¬ 
ues are separated by more than 2 kilometers for 
proton and iron. So, one can use this procedure 
of calculating mean muon production height for 
the primary mass determination. 
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Figure 7. Comparison of the true (zo) and recon¬ 
structed, using parameter £, (MPH) muon pro¬ 
duction height distributions for two primaries, 
proton and iron, of Eq=10 15 eV. Data sample 
was enriched in high energy muons by selecting 
muons with £ < 0 . 02 . 


Parameter £ can be used for investigation of 
muon rapidities in showers. Substituting momen¬ 
tum for the total energy, which is allowed above 
1 GeV without signifcant error, one can express 
rapidity y and pseudorapidity 77 using parameter 
( as follows: 


y = | ln 


E+p || 
E ~P\\ 


1 l n (\/C 2 +l)+l 

2 m (-y/c^+D-i 


rj = I 11 —^ ~ _l n i 

As an example in Fig. |] pseudorapidiy distri¬ 
bution of muons in showers initiated by proton 
primary with the energy Eq=10 16 eV is shown. 


Only muons with momentum larger than 1 GeV/c 
are considered. Dotted line shows all pseudora¬ 
pidity values for all simulated muons. Solid line 
corresponds to the case, where 77 can be calculated 
using parameter £, i.e. for muons with absolute 
values of r and p less than 0.4 radian. Dashed 
line corresponds to the case, where the distance 
of muons to the shower core is in the range of 
20 - 180 m. The large values, being the most 
interesting ones, are cut away due to 20 m min¬ 
imum distance, so it is important to come with 
measurements as close to the core as possible. 



h 

Figure 8 . Muon pseudorapidity distribution. For 
explanations see text. 

The above example shows, that by means of ( 
it is possible to investigate rapidities of muons 
in showers. In general, measuring directional pa¬ 
rameters of muons in showers by means of the </ 
parameter one can investigate muon momentum 
space. 
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